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Abstract 

We compute the parton distributions for the chiral quark model. We present a new tech- 
nique for performing such computations based on Green functions. This approach avoids a 
discretization of the spectrum. It therefore does not need any smoothing procedures. The 
results are similar to those of other groups, however the distributions peak at smaller x. 
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1 Introduction 



Nucleoli structure functions have been measured since three decades, yielding very precise 
and detailed information about the parton distributions. Their evolution with Q 2 is one of 
the important tests of QCD |L] and the predicted asymptotic freedom. There are various 
phenomenogical ansatze for the form of the parton distribution, which necessarily has to be 
chosen at some fixed value of Q 2 . In principle these functions contain important information 
on the structure of the nucleon, which in turn should yield information on the dynamics 
that determines this structure. The problem is that at low Q 2 where one should expect the 
nucleon to be described by quark model wave functions QCD becomes nonperturbative so 
that a description by gluons and current quarks is no longer accessible. Rather one expects 
the system to be described by other effective field degrees of freedom. Among various 
models which are used to describe the low energy properties like form factors, couplings 
to the pion and electromagnetic fields or resonance excitations the chiral quark model 
contains several elements a realistic model should have: there are constituent quarks, there 
are mesons, the model has chiral symmetry. The nucleon is described as a chiral hedgehog 
with an occupied valence quarks bound state, as an extension of the Skyrme model. The 
fact that the model contains quarks, and that the meson field is not elementary, but itself 
a quark-antiquark condensate, implies that the coupling of the electromagnetic current to 
the nucleon is described entirely by the quark degrees of freedom. This property makes it 
a promising issue to compute the nucleon structure functions in this model. 

This has been done recently by several groups [§]-[§] . We do not have anything to add to 
the general approach and its basic formulae relating the parton distribution to the quark 
mode functions or the quark Green function in the external field of the chiral hedgehog. 
We also make a number of technical assumptions which have emerged to be necessary. 
So the pion field is varied only on the chiral circle, the problem of translation variance 
remains, and the Pauli-Villars subtraction with a finite cutoff is applied to the quark sea, 
not to the valence quarks. This is done consistently as well for the self-consistent nucleon 
solutions as also for the parton distributions J/J. 

Where we depart from the previous approaches is in the method used for the practical 
computation of the sea quark distribution functions. All previous approaches start by 
discretizing the mode spectrum, introducing a large "confining" sphere with appropriate 
boundary conditions for the quark wave functions. As the spectrum is discrete, the parton 
distributions are not smooth, and they have to be modified by a smoothing procedure. 
Though this, and the discretization as such, could introduce some arbitrariness the various 
groups obtain similar, though slightly different, results. 

It has been demonstrated previously that such a discretization can be avoided entirely 
when computing the energy of the "quark sea". This quantity has been computed by 
Moussallam || using phase shifts, an approach well-known from kink quantization [[|. 
Alternatively |TI| the self energy can be computed from the Euclidian Green function. 
A self-consistent computation of the hedgehog profile for the nucleon was carried out 
recently | JLl[j . We will present here a technique for computing the parton distribution in a 
continuum approach, i.e., without any discretization. 
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While the self energy can be computed using the Euclidean Green function, by intro- 
ducing and deforming a suitable integration contour in the complex energy plane, we have 
to use the Green function in Minkowski space-time if we want to evaluate the parton dis- 
tribution. So instead of working far from the physical cuts of the Green function, we have 
to work on the cut, here. The technique we use is similar to Moussallam' s approach. But 
of course we compute different physical quantities. So we need more information than just 
the scattering phase shifts. Furthermore we will make use of the techniques well-known 
in potential scattering [12] in order to find expressions for the S matrix elements that 
are suitable both for discussing their analytic properties and regularization, and for nu- 
merical computation. Besides avoiding discretization another advantageous feature of our 
approach is the inclusion of the asymptotic spectrum. Indeed the power behavior of the 
asymptotic tails of the angular momentum summation and of the momentum integrations 
is known and these tails can be included on the basis of fits. 

The plan of the paper is as follows: In section 2 we introduce the model and the 
expression which relates the parton distribution to the quark Green function. A short 
derivation of this relation is given in Appendix A. In section 3 we find analytic expressions 
for the Green function and for the S matrix, for the case of a bosonic quantum field in an 
external field. The corresponding expressions for the fermions, used in the real calculation, 
are given in Appendix B. In Appendix C we derive a relation for the Fourier transform 
of the Green function, which is the basis for the numerical calculation. The numerical 
procedure is described in section 4, in section 5 we discuss the results and present our 
conclusions. 



2 The model 



The Nambu-Jona-Lasinio (NJL) model is defined by the Lagrangian [ 13|1 

G 

£njl = VK*7^ -m)i) + - [{^f + (^75^) 2 



(2.1) 



Introducing Lagrange multiplier fields for the scalar and pseudoscalar densities, the action 
can be written in bosonized form as 



5, 



NJL 



*7^ - g(a + m ■ T75) ip - y (a 2 + tt 2 ) + —^—a J . 



2 m/j, 2 



(2.2) 



The parameters of the model are the fermion self-coupling G, the quark mass m. Fur- 
thermore, as the model is non-renormalizable, a cutoff A has to be introduced. In the 
bosonized version these parameters appear as the quark-meson coupling g and the sym- 
metry breaking mass parameter /1. They are related to the basic parameters as g = /iyG 
and m/i 2 = gf^m 2 . The latter equation expresses /1 in terms of physical constants and of 
the coupling g which remains a free parameter. A further relation is obtained from the 
gradient expansion of the effective action. The resulting kinetic term of the pion field is 
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normalized correctly if the Pauli-Villars cutoff is fixed as 



A = M 



\ 



exp 



' 4tt 2 

nJ 2 



(2.3) 



Here M is the "dynamical quark mass" M = gf n . 

We will not consider the explicit breaking of chiral symmetry by the "current" quark 
mass m, so we work in the limit = 0. Furthermore we will restrict the variation of a 
and 7r to the chiral hypersphere 

a 2 + 7r 2 = ^. (2.4) 

The the second term in Eq. [2.2| becomes a constant that can be chosen to be zero. In 
this CclSG, clS shown by in proper time regularization, and in ]nj in the Pauli-Villars 
regularization used here, the model admits a solution which can be identified as the nucleon. 
The self-consistent nucleon configuration is based on a chiral hedgehog ansatz 



g [a + i^ 5 T ■ 7r] = Mexp {i^T ■ <^(x)} = M(x) 



(2.5) 



with 



0(x) = M(r) (2.6) 

for the a and 7r fields. In this external field the quark field has a bound state solution 
identified with the valence quarks. Furthermore the external field leads to a modification 
of the "sea quark" energy. The self-consistent solution is obtained by minimizing the energy 



I 



£quark = -Tl log -i^dp + M(x 



(2.7) 



where the trace extends over the negative energy levels as well as over the bound state. It is 
understood that the vacuum energy is subtracted. Furthermore a Pauli-Villars subtraction 
using the regulator A is required. 

We use here profiles i?(r) obtained by us previously ||TT|. There the sea quark energy 
was obtained as a trace over Euclidean Green functions. 

It has been noticed by several authors that the computation of the parton distributions 
can be related to the Green function of the quarks in the external field. In contrast to 
the computation of the energy we have to use the Green function at real energy, not the 
Euclidean one. This requires some changes with respect to our previous approach, but we 
will retain its benefits: to avoid working with discrete levels and the corresponding wave 
functions. 

We rederive the expression for the parton distribution in Appendix A, following ||, we 
obtain for the sea quark contribution 

d A k 



Di{x) = -ImV c M N J -—5{M N x-k 3 -k 



xe(-fc - M)tr k(l + 7o7 3 )Mk, k, k ) 



(2- 
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For the antiquarks one has to replace x — > — x and to reverse the sign of the argument on 
the function. Here and in the following denotes the nucleon mass, x is the Bjorken 
scaling variable, and (Sp(k, k, ho) is the Green function of the fermions in the external 
hedgehog field. It is the Fourier transform of the position space Green function 5f(x, y,t) 
which satisfies 

[ny9£ - M(x)] 5 F (x, y, t-t')= 5(x - y)5{t - f ) (2.9) 

with Feynman boundary conditions. As the external hedgehog field is not translationally 
invariant it depends on two separate spatial variables, and so does its Fourier transform 

5 P (k, k', h ) = J dtd 3 xd 3 ye l(kot ~ k - x+k '- y) S F {x,y,t) . (2.10) 

The relation to the positive and negative energy eigenfunctions U a , resp., V a is formally 
given by 

S F (x,y,t) = -i(T(V(x,t)^(y,0)> 

= iQ(t)J2e~ iEat U*(x)U a (y) (2.11) 
«>o 

_, (_ t )^e^lV Q (x)K(y). 

Here a > (a < 0) symbolizes the positive (negative) energy states. We will not use this 
definition for the practical computation, as it would necessitate the discretization of the 
energy spectrum. Rather the Green function is obtained in terms of mode functions as 
described in 

As usual the Green function Sp is written in terms of a bosonic Green function using 
the ansatz 

S F = (E + H)G (2.12) 

where G satisfies 



H 2 - E 2 ) G(x, x', E) = [-A + M + V(x) - E 2 \ G(x, x', E) = -5 3 (x - x') . (2.13) 

Here V is the potential or vertex operator 

V(x) = ?7-VM(x). (2.14) 

G(x, x', E) is the inverse of a symmetric operator, it is obtained by standard techniques 
described in Appendix B. 

In addition to the sea quarks the valence quarks contribute to the structure functions. 
The equation for the K p = + partial wave is 

Vo = wgVo- (2.15) 



-A + V 0+ (x 



The spinor ipo is determined by two radial wave functions h(r) and j(r), and the potential 
V 0+ is a 2 x 2 matrix given in [|TTJ. Defining their Fourier transforms as 

/■OO 

h(k) = / drr 2 jo(kr)h(r) , (2.16) 
Jo 

POO 

j(k) = / drr 2 3l (kr)j(r) (2.17) 
Jo 
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one finds 



. bou , . iV c M N r°° ( 2 2 o M N x-£ bou 



Df m (x) = / dfcfc h(k) + j(kY - 2—^— —h(k)j(k) . (2.18) 

71 J\M N x-E hou \ V fc / 

This expression is equivalent to the bound-state in M. The bound state distribution 
function is convergent. 



3 The Green function 



We present here an explicit expression for the Green function that can be evaluated without 
having recurse to discretization, we also derive some further identities that are useful in 
this context. For the sake of transparency we here consider a scalar field in a spherically 
symmetric background. The much more involved expressions for fermion fields coupled to 
the chiral hedgehog are presented in Appendix B. 

We consider a scalar field with a space dependent mass term m(x) which for |x| — > oo 
tends to M. We furthermore assume that m(x) only depends on |x| = r We decompose 



m 2 (x) = M 2 + V(r) . 
The Green function then obeys the differential equation 

—E 2 - A + V(r) + M 2 ] G(x, x ; ; E) = -<5 3 (x - x') . 
The Fourier transformation of the free Green function is given by 

r rl^h p ih(x-x') 

C(°)(x,x';£;) ' ' ' 



(3.1) 



(3.2) 



(3.3) 



(2tt) 3 k 2 - k 2 + ie ' 

where k 2 = E 2 — M 2 . As is well-know the k integration can be carried out and one obtains 

K e i(K+ie)R 



G(°\x,x';E) 



47T KR 



(3.4) 



with R = yr 2 + r' 2 — 2rr' cos(6 l ). It can be decomposed with respect to Legendre polyno- 
mials using the Gegenbauer expansion. One finds 



GW(x,x';£) = ^^P-Jli^ 



yi{Kr : 



l Ji[ Kr i 



Pi{cos{6)) 



47T 



The spherical Bessel functions are defined as in 



-*«E^r^^(«r<)/i t (1) («'/- )P,{v^{(t)) 



(3.5) 



The imaginary part is given by 



ImG(°'(x,x';£) =-kJ] 



(21 



Air 



-ji(nr)ji(Kr')Pi(cos(6)) . 



(3.6) 



The Green function in the external field is decomposed in an analogous way via 

G(x,x';£) = -z«^2i±ll/ r (K,r < )/+(K,r > )^(cos(e)) . 
■ 4vr 

Here the radial mode functions f l satisfy the differential equation 
/ d 2 



2d 1(1 + 1) 



K 2 \f l (K,r) = -V(r)f l (K 1 r) 



\ dr 2 r dr 

The boundary conditions are, in analogy to the spherical Bessel functions: 



limf (K,r)/ji(Kr) 



lim f + (K,,r)/h\ 1 \nr) 



const. 



1 . 



(3.7) 



(3.* 



(3.9) 



The normalization of f l (n, r) is determined by the Wronskian 



(3.10) 



It is useful for the numerical computation as well as for the separation (subtraction) of the 
free Green function to split off the Bessel functions via 



/f («>r) 



l + cj(K,r) 
1 + c[{K,r) 



h 



(i), 



nr) 



ji(Kr) 



(3.11) 
(3.12) 



The second equation is used only below the first zero of ji(nr). In fact the function 
c{(k, r) is pathological as it has poles and zeros at the zeros of ji and ff, respectively. An 
alternative, more useful expression for ff is obtained by forming a linear combination of 
f^f and (/z + )* which becomes regular as r — > 0; explicitly 



ff(K,r) = U(l + 4(K,r))h\ 2 \ K r) 

+ irJS( 1+ ^-))^w}- 



(3.13) 



The function (k, r) tends to zero as r — > oo and goes to a constant as r — > 0. It is well 
suited for numerical computation, using the differential equation obtained by inserting the 
ansatz (|37TT|) into Eq. ( ggfl . 

The function cf is related in a simple way to the S matrix. The relation follows from 
considering the regular solution ff in the form ( |3.13| ). Indeed the S matrix is defined via 
the asymptotic behavior of the regular solution: 



fl 



hf\nr) + e 2i5 ^ E) hf\^r) 



(3.14) 



h 



<~W^ , 1 + ^( K > 0)^(1), 



[Kr\ 



l + c+(«,0) 



nr) 
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So the S matrix is given by 



Sl{E) ~ e -l + c+(«,0)- (3 - 15) 



As is well-known ||, the phase shift can be used to compute the zero point energy. For 
the quark sea in the chiral quark model this has been done in ||; the relation is, in this 

case, 

Eo = -- V(2^ + 1) r dnVn 2 + M 2< ^4^ . (3.16) 
7r K p Jo an 

For the parton distribution we need the imaginary part of the Green function in mo- 
mentum space at equal three-momentum. We have 

ImG(k,k,£) = J <?x J dVe" ik(x - x,) i [G(x, x', E) - G(x', x, E)] . (3.17) 

As the Green function is a symmetric in its spatial arguments the expression in brackets is 
just twice the imaginary part of the Green function in position space. It is straightforward 
to show that this imaginary part is given by 

Im G(x, x', E) = -kJ2 («, r)fr («, r)P,(cos 6) . (3. 18) 

i Alx 



This is analogous to the imaginary part of the free Green function ( p.6|) . Moreover the 
Fourier transform factors in a similar way. If decomposed into partial waves the Fourier 
transform reduces to a Fourier-Bessel transform and one obtains 

Im G(k, k, E) = -4™]T(2Z + l)/f («, k)fr(K, k) (3.19) 

i 

With ,00 

f l -(n,k)= drr 2 ff(K,r) 3l (kr) . (3.20) 



This Fourier transform is problematic as the integral is not absolutely convergent. This 
problem is analyzed in Appendix C, making use of Wronskian identities for the free and 
exact Green functions. The final result is 

/,"(«,*) = (e^ + l)^-{(-l) l+1 5(k + K)+5(k-K)} 



+ ^71F / drr a j z (Ar)V(r)/f( K ,r) . (3.21) 



Akn 

, r 

K 2 — k 2 Jo 

Within the range of integration for the parton distribution, Eq. fl2.8j) , the 5 functions do not 
contribute as k 2 — k 2 > 0. So in fact only the last term contributes, the integral converges 
as V(r) — > exponentially as r — > oo. This provides a numerically stable expression for 
this Fourier transform. 
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4 Numerical computation 



In the main text we have considered scalar fields in order to present the scheme for evaluat- 
ing the parton distributions on the basis of Green function techniques. We will continue to 
do so here. Obviously the real computations were done using the formalism based on the 
Grand spin reduction of the fermion Green function as presented in Appendix B, leading 
to 4 x 4 complex coupled differential equations. 

The functions cI(k,t) have been calculated by solving the differential equation 

d 2 (l h\ lY (Kr) 



dr 2 




( K ,r) = -V(r) 1 + 4(K,r) . (4.1) 



which follows from Eq. ( |3.8| ) by inserting the ansatz ( |3.11| ). We have used a simple four- 
step Runge-Kutta-scheme, the integration was started at values of R much larger than the 
range of the potential, where c z + vanishes by definition, i.e., by its boundary condition. The 
accuracy of these solutions was checked by the Wronskian relation, which was constant to 
at least six significant digits. As a further check we used the unitarity of the S matrix as 
the 4x4 matrix relation (p,13|) . This is displayed in Fig. [I] 

The functions ff were then obtained in the form (|3.13| ). As the Hankel functions 
become singular at r = this form is not suitable at very small r. In this region we have 
used the form (|3.12 ), starting the integration of the differential equation at r = 0. This 
solution was normalized correctly by fitting it to the composed solution at some value or 
r where both solutions are reliable. This is displayed (for one component of the real 4x4 
solution) in Fig. |^. 

The next step in the computation of ff{n, k), i.e., the Fourier-Bessel transform, using 
Eq. ( |3.21| ). The Bessel functions were generated recursively, for small r they were ob- 
tained via power series expansion. The integrand is then obtained by squaring the Fourier 
transform and by supplying the pole prefactor l/(/t 2 — k 2 ) 2 . In the fermion case the cor- 
responding operation is given in Eq. ( |B.14j ). In this way the integrand is evaluated on 
a lattice in the variables k and k. The region of integration is displayed in Fig. |3|; the 
dependence on the scaling variable x comes from the boundary k > k min = \xM^ — fc | 
with k = —\/ k 2 + M 2 and some simple prefactors, see Eq. ( B.19|) . The first step is the in- 
tegration over k. As we work with a discrete lattice the lower end of integration could lead 
to a step structure; this is avoided by fitting the behavior near k min down to this boundary. 
The integration over large k is unproblematic, both due to the prefactor 1 / (k 2 — k 2 ) 2 as to 
the decrease of the squared Fourier transform. 

The next step is the k integration. Here the lower boundary is k = and the integrand is 
well-behaved there. The behavior at large k is powerlike, as expected in a loop integration. 
The integration over the computed integrand was extended to values k ~ (25 -f- 35)M, 
depending on K. The integrand decreases as k~ a , so the integration up to oo is appended 
using an asymptotic tail of the form Ak~ a + Bk~ 5 . An example is displayed in Fig. |]. This 
integration over k has to be performed at each value of x. 

The final step is the summation over Grand spin K and over parity P. The behavior 
of the terms in this sum is displayed in Fig. H. We plot the terms for the different parities 
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for the physical quark mass M and for the cutoff mass A. The K" 1 behavior reflects the 
divergence to be expected from second order perturbation theory. The difference of the 
quark mass and cutoff is taken in the form J7| 



"V, 



A 2 



x,A) 



(4.2) 



These terms decrease as K~ 3 . The sum over the K-spin was extended to K — 18. The sum 
over the higher angular momenta was appended using a power fit to the terms computed 
numerically. 

We plot the final results for the isoscalar unpolarized distribution in Figs. |7] and |[ The 
antiquark distribution is obtained as Di{x) = —Di(—x) 

It is interesting to compare these results with second order perturbative contribution 



DAx) 



d 3 k 



x 



4tt 2 J (2tt) 3 
A 2 + k 2 



6{k 3 



\x 



M N )Tr (u(k)[U(k)\ 



In 



k 2 (A 2 -M 2 



with 



M 2 + K 2 (A 2 + K 2 )(M 2 + K 2 ) 



\x\M N (k 3 - |x|M N )k 2 



(At 



3^2 



(4.3) 



(4.4) 



This is done for the sum of quark and antiquark distributions in Fig. |6|. It is a further 
numerical cross check to compute the second order perturbative contribution form the 
mode expansion. This can be done by replacing in the expression for the Fourier transform 
of the mode function f^~{n, k), Eq. ( |B.15| ), the exact solution r) by the free solution 

^nJK n (^r). This comparison is included as well in Fig.||. The agreement is not perfect, but 
it proves the overall consistency and illustrates the precision of our results. As the higher 
order perturbative contributions, and therefore their error, are less important this graph 
can be used to improve the overall result by replacing the computed second order by the 
analytic one. 

Finally, it was also checked that the relation 



f L (A(z) + Di(-x))dx = N c 
Jo 



(4.5) 



was fulfilled within 0.1%. In principle this relation || is valid exactly if the integration is 
extended to x = oo. As seen in the figures, the integrand is already very small at x — 1. 



5 Results and Conclusions 

We have presented here a new evaluation of the isoscalar nucleon structure functions, using 
the continuum approach based on the quark Green function in an external background 
field. This constitutes a natural extension of our use of Green function methods to the 
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computation of self-consistent solutions ||11||. The method presents various internal cross 



checks as mentioned in the previous section. It is a relatively fast method, the computation 
of the structure function takes about 5 hours on a standard PC with a 450 MHz Pentium 
processor. 

The self-consistent solution was computed using our previous methods. For the quark- 
meson coupling we have chosen the value g — 4 which was found to be preferred by the 
comparison of the predicted static parameters with experiment. For f n we have taken the 
experimental value 93 MeV. 

The results for the parton distributions are shown in Figs. ^ and || As compared to 
previous analyses, based on the explicit wave functions of a discretized system, we find an 
overall agreement as far as order of magnitude and typical structure are concerned. The 
results are similar to those of Weiss and Goecke and of Wakamatsu and Kubota 
There is a tendency, however, of our structure functions to peak at somewhat smaller values 
of x. In view of some differences in the self-consistent hedgehog profiles and parameters 
this is however a very satisfactory agreement of entirely different numerical approaches. 

Our results are compared to the parton distributions obtained in the GRV parametri- 
sation [|Tj]] at Q 2 = 0.40 GeV 2 , that is, of the order of the cutoff A 2 . The agreement is 
certainly qualitative and even semi-quantitative. 

One of our main results that our method corroborates previous analyses in a very 
satisfactory way. However, it also adds a new and efficent tool for such and similar com- 
putations. As we have shown it comprises various internal cross checks, some of which 
are due the use of analytical methods of potential scattering. Furthermore the method is 
comparatively fast, and so can be useful for more elaborate computations. 

The most obvious continuation of the present work would be the computation of the 
isovector structure functions. As it is well-known these functions vanishe in order , and 
the first nonvanishing part is of the order ~ fl 1 . These parts are necessary in order 
to obtain the u- and d- quark distributions separately. In analogy to the computation of the 
momentum of inertia |TT| (first done by Wakamatsu and Yoshiki |18| using discretization), 



it is feasible to extend our method to A^-distribution functions, taking into account the 
rotation of the soliton. 

The fundamental formulae have been found after a long discussion in ||19|| , || and p0 |. 



Using the same procedure as for the iV^-case, they can again be written in terms of 
Minkowskian Green functions. One obtains 

W*> -'<*)]<" = - 2 ^I.n/|^/(g;*o + ^-.M N ) 

xTYS F (A, k', A„)T a 7 S F (<:', k, k^rj^' + 7 3 )e(-fc°) (5.1) 

and 

[u(x)~d(x)f ) = - 2 ^^-^Mx) + d(x)\, (5.2) 

The computation of the fermion Green functions in the external hedgehog field has been 
outlined here. Comparing to the computation of the isosinglet structure functions we need 
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one more integration, the one over k'. In order to avoid Fourier transforming it can be done 
in x space. Furthermore one looses the factorization property (|3.19 ). On the other hand 
the convergence in momentum and angular momentum will be faster, as for the analogous 
leading order Feynman graph. Furthermore the functions themselves are smaller, so one 
can somewhat reduce the numerical accuracy. So this computation seems feasible and is 
being undertaken at present. 

A further issue will be the computation of the polarized structure functions 0,[5| and 
of the skewed parton distributions [21|. It would be interesting, furthermore, to study the 
full Q 2 dependence, as has been done for the bound state contribution by Weigel, Gamberg 



and Reinhard p2| as well as by Ruiz Arriola [23 



As a side result we have obtained a very fast method for calculating the phase shift. 
Alternatively to the method employed in P,|24]|, the computation of the system of smooth 



nonoscillating functions c" + allows the extraction of the multichannel S matrix. A further 
possibility - and cross check - is the use of (CTTB) and its multichannel extension. 
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A Parton distribution 



The fermion Green function has been formally expressed by the eigenmodes U a and V a in 
Eq. ( [2. 11| ). Its Fourier transform ( |2.10| ) is given by 



U a (k)U a (k>) 



V a (k)V a (k') 



swk, k ; , kn) = — — — , 

where the Fourier transform for the eigenspinors is defined via 

d 3 k 



UJx) 



(2tt) 3 
d 3 k 
(2tt) 3 



V a (k)e 



ikx 



U a {k)e 



ikx 



The parton distribution of sea quarks is given by 



q[x) 



iV c M N 
2tt 



(A.l) 



(A.2) 
(A.3) 



<f'X I dz°e ixM " z ° Y, e i|£Q|z Vi(-X)(l + 7 o7 3 )K,(-X - z°n 3 ) . (A.4) 

00 a<0 
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If one rewrites this in terms of the Fourier transforms, the z° and X integration can be 
performed with the result 

d 3 k 



q(x) 



N r M] 



N 



(2tt) 3 



Y,S(\E a \ + xMx-k 3 )tr 

a<0 



(l + 7o7 3 )K(k)^(k) 



(A.5) 



We note that 1 + 707 3 is a hermitean operator. Considering the Green function in the 
form (|A.1|) we see that 

^ F (k,k,A; )7 -7 4(k,k,A ;o ) 

= £ 2nS(k° - E a )U a (k)Ut(k) - £ 2*5(k° + \E a \)V a (k)V^(k) . (A.6) 



Q>0 



o<0 



We therefore find 



q(x) 



-Im N r M] 



N 



d 4 k 

W) 4 



5{xM N - k 3 - k )Q(-k°)tr [( 7o + 7 3 )5 F (k, k, A; )] . (A.7) 



The restriction of the integration to negative k° does not appear in (A. 9) of 0. It does 
appear in their Eq. (7.3) for the second order perturbative contribution. As is apparent 
here it is crucial for restricting the integration to the sea quark states. We have remarked 
below Eq. ( |3.21|) that k 2 — k 2 > so that the 8 functions do not contribute. In terms of 
the variables used here 



k — k = k — (k°Y + M > (k ) + M — (ky = M + x 2 M^ - 2xM N A; u > . (A.8) 

The antiquark distribution is determined formally by the positive energy eigenspinors. 
This means that one has to replace the G function by Q(k ). At the same time the energy 
of these levels has to be replaced by —E a , a > 0. In the Green function this reversal of 
sign is obtained by replacing x by —x. 



B The Green function for the quark system 
B.l The boson Green function 

As we have described in section 2 the fermion Green function has been reduced to a bosonic 
one which satisfies the differential equation (2.13). It is a 8 x 8 matrix. Using parity it 
becomes block diagonal with two 4x4 matrices. As the system is invariant with respect 
to K spin, it can be expanded [^5| with respect to K-spin harmonics ^n' Kz y i a 

G(z,z',E)= £ g«£(r,r', K )E«> K *(z)®E^(z') (B.l) 

K,K Z ,P 

with k = \/E 2 — M 2 . The radial Green functions are 4x4 matrices. They can be written 
in terms of 4 - component mode functions as 



T.T , K) = K 



e(r-r')f« + (K,r)fr(K,r') (B.2) 



+0(r'-r)f«-(K,r)fr(K,r') 
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The mode functions satisfy the differential equation 



d 2 _2d_ K n (K n + 1) 



dr 2 r dr 



(B.3) 



The subscript denotes the four components, the superscript a = 1 ... 4 labels the four 
independent solutions which form a fundamental system. The potential V™ (r) is the grand 



spin partial wave reduction of V(x), Eq. ( 2.14 ). Its explicit form has been given, e.g., in 
Appendix A of [|TTJ . 

In order to split off the behavior at r = and at r = oo of the free solutions we write 
them in the form 



f 



(B.4) 
(B.5) 



Here K m takes the values K — 1, K, K, K + 1 for the four components labelled by m. As 
for the single component case presented in section 3 the boundary conditions are 



lim c° + (k,t) 
lim c£~(K,r) 



0. 

const. , 



(B.6) 
(B.7) 



and the fundamental system /" is fully specified by the Wronskian 



r 2 W^(f + ,f- 



J n 



dr 



f~-f 

J m J i 



d 



ra - 

"'■ dr Jm 



■_ 5 a0 



(B.8) 



As for the single- component case /° can be written as 



fi'M = i|(^ + ^ + («,0)) [(5 + c + («,0))~ 1 ]"(^ + ^+(«:,r))/ i «(«:r) 

+ (tfi + 4 + («.*-))fcg,(«r)} (B.9) 



and it was computed in this form except at very small r. The 5 matrix is defined by the 
asymptotic behavior of the regular solution r), from the representation (|B9|) we 

find 



= (i + c + M)) ([i + c+M) 



(B.10) 



From this form it is not obvious that it is unitary. We have to use the further information 
that the Green function forms a symmetric matrix 



g m n(R, R, «) 



9nm(Ri R, 



(B.ll) 
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which leads to the relation 

(1 + c(k, 0)) n a [(1 + c(k, 0))- 1 ] a = (1 + c(k, 0))™ [(1 + c(k, 0))- 1 
Using this relation one obtains 



(B.12) 



S^S* x n = (i + c+(k,0) 



1 + c+(k,0) 



l+C+(/€,0))^[(l+C + (/€,0)) 1 

(l + c + («, 0))* [(l + c + (k, 0)) _1 j ° (l + c+(«, 0))" [(l + c + (/c, 0)) _1 



= 5 m „. (B.13) 

This relation can be used to check the numerics. The inaccuracy is less than 10~ 6 . Using 
the symmetry relation (|B.11|) in the form (B.12| ) one can show that the imaginary part of 
the Green function again factorizes: 



lmg mn (r,r',K) = -Kf° (K,r)/£ («,r'). 



(B.14) 



For the numerical computation one needs the Fourier transform of the mode functions 

fm~- lt is s iven b y 

fZT{K,k) = (e 2i5 ^ + l) a J-6(k-K){(-lT +1 6(k + K)+8(k-K)} 
V 



K 2 — k 2 Jo 



drr j Km (kr)V mn {r)f* («,r) 



where e 2lS ^ + 1 means 



(e 2i5 W + 1 



(B.15) 
(B.16) 



With these preparations we can obtain the Green function for the fermion system. We 
use the Fourier transform of Eq. ( |2.13| ): 



5(k, k, k) 



-Vk 2 + M 2 - kak + 7o M) G(k, k, 



The action of ctk onto the spinor harmonics is given by 



erk 



( H!(*) ^ 






Sa(*) 




Si(*0 


3s(*) 




24(fc) 









(B.17) 



(B.18) 



Combining the H spinors, the 7 matrices and Eq. ( P . 1 7|) one finds after some algebra 
|dQ k Imtr [(70 + 73) £(k,k,K)] = Im £(2if + 1) j (2£ - zM N + m) (# u + £44) 



2£ - xM N - M ) ( #22 + (733 



(B.19) 



-2(xM N - E)-[g 12 + g 3 4 ) + 2fc( #i 2 + £34 ) • 
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Here we have suppressed labels and arguments of gij; explicitly they are given by 

1^ POO 

lmg^ P (k,k,K) = -— —2 / drr 2 j Ki {kr)V in {r)f£-(K,r) 

(k 2 — k 2 ) Jo 

00 

2 



x/ drr 2 j Kj (kr)V jn ,(r)f*r(K,r). (B.20) 



C Fourier transform 

We here discuss the Fourier transform of the mode function ff(n,r) 

poo 

fC(K,k)= / drr 2 ff(K,r) 3l (kr). (C.l) 





This Fourier transform is problematic numerically as the integral is not absolutely conver- 
gent. 

We start with the differential equations for ff(n,r) given in Eq. ( |3.8|) and for ji(kr) 
which reads 

'^r^j l (kr) + l ^-j l (kr) = k 2 j t (kr) . (C.2) 
r z dr dr r z 

From these we derive 

- J o dr ^K^)^ r2 ^/f («> r ) - ^~( re > r )^ r2 ^(^ r )J 
f R ( \ f R 

+ J drr 2 ji(kr)V(r)ff(K,r) — yK 2 — k 2 J J drr 2 ji(kr)ff(K,r) (C.3) 
and after integration by parts 

rR l ( rR 

drr 2 ji(kr)ff(n,r) = — —I / drr 2 ji(kr)V(r)ff(K,r) (C.4) 

fir — K Jo 



-R 2 



ji(kR)±fr(K,R) - fr(^R)^UkR) 



The boundary contribution at r = vanishes, because both ff and ji behave as r . The 
first term on the right hand side is well-defined in the limit R — > 00, due to the decrease of 
the potential. The two terms in the bracket display an oscillatory behavior with constant 
amplitude due to the asymptotic behavior of ff : 

/f for) * /f°°(«,r) = \ (fcf^rj + ^^r)) . (C.5) 

This behavior makes the limit R — > 00 of the integral on the left hand side ill-defined. It 
is these asymptotic oscillations that we have to split off. 
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Using the Wronski determinants for the free Bessel functions and Eq. ( |C.4j ) one finds 



R 



drr 2 ji(kr) (/, (n,r) - f l 
1 



Ik, r 



k 2 - k 2 



R 



drr ji(kr)V(r)f l (K,r) 



-R 2 



d 



ji{kR)—{fr{^R)-fr°°{^R) 



-(fr(K : R)-fr°°(K,R)) -^Ji(kR) 



(e 2 ^ - l) 



(C.6) 



The limit R — > oo of this equation is well-defined. One obtains 

1 



drr 2 ji(kr) (j t (k, r) - f l '°°(K,r)) 



k 2 — k 2 
k l 

—i 



drr ji{kr)V{r)f l (n,r) 



£ 2iS K _ l 



(C.7) 



Both parts of the integral on the left hand side are still ill-defined. 

In scattering theory this problem is solved by working with normalizable wave func- 
tions obtained by giving small imaginary parts to the momenta (here k) and taking the 
appropriate limit to the real axis. This is of course the origin of the physical cuts. This 
is usually done using Jost functions flI5| , like / ; + (/t,r), which asymptotically behave as 
exp(±mr). However, for brevity of presentation, we have preferred to work here with ff 
directly. The following results are a heuristic presentation of what we have obtained by 
using the meticulous ie prescriptions. 

ji(nr) + iyi(nr) are 



The Fourier transforms of the free solutions h\ (kt) 



drr z 



MkrMKr) = ^- {(-l) l+1 5(k + k) + 5{k - «)} . 



and 



drr 2 ji{kr)yi{Kr) 



V 



[k 2 - k 2 )K l + l 



(C.9) 



One obtains 



fr(K,k) = (e 2 ^ + l)^-{{-l) l+1 5{k + K)+5{k-K)} 
V 



+ 



K 



k 2 Jo 



drr 2 ji(kr)V(r)f l (k, r) . 



As we have mentioned in the main text, the distribution character of this result does not 
play a role in the numerical computation. The finite integral which appears on the right 



means principle value 
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I 

K JO 



hand side plays a central role in evaluating the parton distributions. Using further the 
methods of potential scattering we can derive a useful identity relating this integral to the 
S matrix. 

The integral equation accompanying the differential equation (|3.8|) for the solutions f + 
reads 

fi~(n,r) = —in I dr'r 2 ji(Kr')h\ 1 \nr)V(r')f l h (K,,r') 
Jo 

-in j°° dr'r'^j^Kr^iKr^Vir^f+i^r') . (C.IO) 

Jr 

As asymptotically ff(n, r) behaves as h^far) one finds 

dr'r' 2 ji(Kr')V{r')ff{K,r') . (C.ll) 
Combining / _ using / + in the usual form 

/f(M = l(rt(^r)+e^f+(K,r)) (C.12) 
the integral takes for k = k the form 

/ drr 2 j l (Kr)V(r)fr(K,r) = — (e 2 ^ - l) . (C.13) 

This relation has been used as a useful crosscheck. A further cross check is obtained form 
the Born approximation to this relation: 

- k r drr 2 jf(Kr)V(r) = (C.14) 
Jo 

which determines the behavior of the phase shifts at large I and at large k. 
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Figure captions 



Figure 1: Unitarity of S matrix: (S'S) \(k) for K p = 1 + : solid line: total matrix element, 
other lines: the eight components Sl a Si, a = 1 . . .4, product of real parts and product of 
imaginary parts. 

Figure 2: Behavior of fi~(n,r) for (K p = 5 + , k = 1), at small r. dashed line: computed 
via Eq. (|3.13| ), long dashed line: computed via Eq. ( |3.12|) . 



Figure 3: Integration area in a k — k plane: the dashed line represents position of the 
pole, the lower solid line is the minimal lower limit of the integration, the upper solid line 
shows the lower limit for a positive value of x. 

Figure 4: Convergence of the K integration: integrand for the mode function at x — .25 
for K p = 5 + ; solid line: the integrand, dashed line: power fit with Ak~ 4 + Bk~ 5 . 

Figure 5: Behavior of Grand spin contributions to Di(x = .25) at large values of K: solid 
and dashed lines: parities + and — with dynamical mass; dotted and long-dashed lines: 
parities + and — with cutoff mass; dash-dotted lines: regulated contribution. 

Figure 6: Second order contribution; solid line: computed via a Feynman graphs; dashed 
line: computed via second order mode sum. 

Figure 7: The isosinglet unpolarized distribution of quarks and antiquarks | [q{x) + q(x)] 
solid line: total result, dotted line: valence contribution, dashed line: sea contribution, 
squares: NLO GRV [0] parametrisation 

Figure 8: The isosinglet unpolarized distribution of quarks and antiquarks ~ [q(x) — q(x)] 
solid line: total result, dashed line: sea contribution, squares: NLO GRV |~7] parametri- 
sation 
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